USRARY 

ns naval POSTGR^DU ATE SCHOOl 

MONTIRL'^. CALlfOmiA 



r 



i: 







# 



r 








T 



I 




■’1 



AN INVESTIGATION OF A METHOD FOR 
UNIQUE DETERMINATION OF AIRCRAFT 
LONGITUDINAL STABILITY DERIVATIVES 
FROM TRANSIENT FLIGHT TEST DATA 



R. R. Cornwell. Lt. USN 
F, F. Eckhart 

Aeronautical Engineering Report No. 459 



Submitted in partial fulfillment of the requirements 
for the degree of Master of Science in Engineering 
from Princeton University, 1959. 



ACKNO W LEDGEMENTS 



The authors wish to exprees their appreciation to Professor 
Edward Seckel for his considerable assistance throughout the conduct 
of this investigation. 

Appreciation is also expressed to the Aeronautical Engineer- 
ing Department of Princeton University for the use of their aircraft and 
flight test facilities. In particular, gratitude is expressed to Mr. Robert 
F. Cooper for his cooperation in maintaining the aircraft and to Mr. J. 
Richard Whitley for the installation of the aircraft instrumentation and 
for his assistance in operating the telemetering ground station. 

The second author wishes to express his appreciation to the 
Institute of the Aeronautical Sciences for supporting his study at Princeton 
University under their Flight Test Engineering Fellowship program. 



TABLE OF CONTENTS 



Page 

Summary 1 

Table of Symbols 2 

Introduction 4 

Equipment 6 

Adr borne Instrumentation 6 

Ground Station Equipment 9 

Instrumentation Calibration 10 

Procedure * 14 

Discussion and Results 16 

Theory 16 

Linear Dependence 21 

Preliminary Investigation 23 

Amalog of the Test /Jrcraft 24 

Solution for Stability Derivatives from Analog Data 24 

Results of tlie Preliminary Investigation 27 

Results of the Flight Test Data 28 

Proposed Improvements 32 

Conclusions 36 

Recommendations 37 

References 38 

Tables 39 

Figures 



44 



SUMMARY 



Using a North American Navion as a Qight test vehicle, a study 
was made of the practicality of experimentally determining all longitu- 
dinal stability derivatives uniquely by releasing a large weight in flight. 
The solutions for the stability derivatives were attempted by matching 
the transient response of the aircraft with solutions to the equations of 
motion generated on analog computer. 

Due to instrumentation difficulties definite verification of the 
practicality of this procedure was not achieved. However, the experi- 
mental results indicated that the stability derivatives could be found 
uniquely by releasing a large weight in flight. 

Unique longitudinal stability derivatives were found from simu- 
lated data generated on the analog computer. Better accuracy for these 
derivatives was achieved if the transient response to an elevator step 
was analysed in addition to the weight release response. 

It is considered that this procedure for the unique determination 
of longitudinal stability derivatives is suitable for further investigation. 
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INTRODUCTION 



In recent years, dynamic flight testing has received progres- 
sively increasing emphasis. The objectives of dynamic testing are 
several In number. However, receiving the greatest attention is the 
determination of stability derivatives from the transient response of an 
aircraft to an applied forcing function. There are many advantages to 
this type of flight testing. Among the more important advantages are 
the short flight time and relatively simple flight techniques required. 

In many cases, dynamic flight testing offers the only means of finding 
certain stability derivatives. 

The basic approach in using dynamic flight testing to determine 
stability derivatives is simply described. A time history of the air- 
craft's transient response to some input forcing function is recorded. 

It is assumed that the form of a set of linear, constant coefficient, dif- 
ferential equations describing the system is known. The time history 
is then used to evaluate the constant coefficients in these equations. 

A basic problem that is present in dynamic flight testing is the 
unique determination of all longitudinal stability derivatives. The most 
common procedure has been to assume some theoretical relation be- 
tween these derivatives. This report is concerned with the unique de- 
termination of all longitudinal stability derivatives without any such 
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asBumptioDB. This is accomplished by using a weight release as a 
forcing function to produce the aircraft's transient response. This 
relieves any linear dependency between the various longitudinal stability 
derivatives. 

This investigation was carried out at the James Forrestal 
Research Center, Princeton University, Princeton. New Jersey, during 
the 1958-1959 academic year. 
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EQUIPMENT 



Airborne Instrumentation 

The test aircraft employed was a North American Navion, 
which is a four place, all metal, low wing aircraft. The Navion is 
powered by a Continental E-185 engine rated at 185 HP at 2300 RPM 
and 29 inches manifold pressure. A three view drawing of the aircraft 
is shown in Fig. 1, and the general specifications are recorded in 
Table I. 

The aircraft was modified by the installation of a Navy Mark 
51 bomb rack which was placed on the lower center line of the fuselage, 
aft of the nose gear wheel well. The bomb rack was installed so that 
the release mechanism could be actuated either electrically by a sole- 
noid or manually by a cable. A 566. 5 pound solid steel bar was equipped 
to hang from the bomb rack. This steel bar was 9. 5 inches in diameter 
and 27 inches long. Henceforth in this report, the steel bar will be re- 
ferred to as the bomb. Figs. 2 and 3 show the bomb installed on the 
aircraft. 



Standard aircraft instruments were used to determine altitude 
and airspeed. In addition, the aircraft was instrumented to measure 
normal acceleration, pitch rate and angle of attack. These variables 
were fed to a telemetering system and recorded at a ground station. 
Also, the instant of bomb release was recorded. Fig. 6 is a schematic 
wiring diagram for the instrumentation system. 
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The accelerometer used was a Giannlni Model 25113-1-20 with 
a range of 1 g<. The sensitive axis of this instrument was canted 
laterally at 45°. With this arrangement, the accelerometer will only 
respond to a 45° component of normal acceleration. This effectively 
multiplies the range of the instrument by the square root of two for 
normal accelerations. The full scade range is therefore - 1.41 g's. 

An initial installation with a > 3 g's Genesco accelerometer was made, 
but unsatisfactory data was obtained. No other small range accelerome- 
ter was available which could be employed with the telemetering equip- 
ment. A i 1 g accelerometer canted lateredly at 45° will of course 
respond to any lateral accelerations. However, the instrument provided 
reasonably good results. Full scale voltage for the accelerometer poten- 
tiometer was provided from the 28 volt battery bus. The resistance of 
the potentiometer is 2000 ohms. A bias battery was included to adjust 
the steady state voltage to the desired level. 

The pitch rate was measured with a Honeywell gyro Air Force 
Part No. JG 7005 A-24. This type instrument operates with a rotor 
power supply of single phase, 400 cycle, 115 volt alternating current. 

The transducer is of the potentiometer type with a resistance of 530 
ohms. The axis of the gyro was carefully checked to align with the air- 
craft axis. The potentiometer was supplied with five volts from the 
battery bus. These same five volts were supplied to the telemeter equip- 



ment as a full scale reference. 
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The angle of attack vane was mounted on a four foot boom which 
was located on the left wing tip. A balsa vane with a direct connection to 
a helipotentiometer was employed to reduce the moment of inertia and 
friction of the system. The potentiometer was supplied with 28 volts and 
a bias battery. 

The telemeter unit was the pulse width type and was built by the 
Applied Science Corporation of Princeton. This unit contains 43 channels 
and a normal sampling rate of 20 per second on each channel. Three 
channels were used for angle of attack and pitch rate while six channels 
were used for the normal acceleration. 

The telemetering unit is designed to use five volts as a full scale 
reference voltage. When properly instrumented, these same five volts 
should be the full scale reference for each transducer. Small variation of 
the aircraft voltage will not effect instrument calibrations with this 
arrangement. However, it was desired to improve the sensitivity of the 
angle of attack and normal acceleration measurements; therefore. 28 
volts were employed as full scale. Since the telemeter will accept only 
zero to five volt signals, a bias battery was employed to set the steady 
state potentiometer inputs within this range. This has the disadvantage 
that a change in the 28 volt bus will change the steady state output of the 
angle of attack vane and the normal accelerometer. However, the sensi- 
tivity will not be effected. 
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Ground Station Equipment 

The ground station equipment included a telemeter receiver, a 
magnetic tape recorder, ein analog computer and a Sanborn model record- 
er. These components are shown in Figs. 4 and 5. 

The telemeter receiver was built by the Applied Science Corpor- 
ation of Princeton. This unit is designed to receive the coded output of 
the telemeter transmitter, decode the signals, mix channels as desired 
and transmit the data to a recording device. 

An Ampex Model 309C dual track tape recorder was employed to 
record the data. This provided a convenient means of retaining data until 
ready for further analysis. The fact that the data was in an electrical 
form was a great advantage when the use of an analog computer is contem- 
plated. All voice transmissions were also recorded on the second tape 
channel. 



A Goodyear Aircraft Corporation Model L 3 (GEDA) linear 
electronic differential analyser was employed to anedyse the data. This 
analog computer provides twenty-four stabilised DC amplifiers with an 
open loop gain of 5 x 10 . The computer had a guaranteed accuracy of 
one percent. Provisions for accurately setting the computer board poten- 
tiometers was available. 



A Sanborn Model 154-100 B four channel recorder was used 
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for visual presentation of the data. This instrument had a high natural 
frequency, and the low frequency response was flat to zero. 

Instrument Calibration 

Three components of the instrumentation system required call* 
bration. These were the normal accelerometer, pitch rate gyro and angle 
of attack vane. The eiircraft airspeed Indication emd sdtimeter were not 
calibrated as previous investigation indicated small errors in these 
systems. 

The accelerometer was calibrated by tilting it through 90®. 
The angle of tilt and the percent of full scale voltage v/ere recorded. 
This checked very closely with the factory calibration which was there- 
fore accepted as valid. Fig. 7 shows the sine of the angle of tilt versus 
the potentiometer output in percent full scale voltage. The slope of this 
curve was found to be 2.113 g's for full scale voltage. Since the instru- 
ment was canted at 45®, this value was divided by the sine of 45® to 
give 2.995 g’s for full scale. As installed in the aircraft 28 volts were 
applied across the accelerometer potentiometer with a five volt portion 
of this representing the telemetering full scede. As received at the tele- 
metering grovtnd station, full scale was therefore: 5/28 x 2.995 > 

. 534 g's. 

Part way through the flight test program, it was discovered 
that when releasing the bomb electrically, the bomb release solenoid 
reduced the battery bus voltage by .6 volts. This reduced voltage lasted 
as long as the spring loaded bomb release switch was held down. This 
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effectively reduced the accelerometer iiero by 12% of telemetering full 
scale. This is indicated by the negative peak on the flight 4 accelerom* 
eter trace in Fig. 14. 

Henceforth, the bomb was released manually. However, in the 
two additional drops which were accomplished the accelerometer was 
found to have a dead zone in the steady state position. Time would not 
permit additional drops. Careful recalibration of the instrument showed 
it was linear outside the dead zone, and the zero was displaced . 6 g. 

The new calibration is shown on Fig. 8 as .518 g'e for telemeter full 
scale. 



The dynamic characterisitics of the accelerometer were deter- 
mined as follows. The instrument was placed on a board with the sensi- 
tive axis vertical. The support was then sharply removed, and the tran- 
sient response of the instrument recorded during its free fall. A drop of 
approximately four inches was sufficient for the instrument to reach a 
zero g steady state. The instrument dynamic characteristics were de- 
termined by matching an analog computer trace to the instrument re- 
sponse. The damping ratio of . 481 and the natural frequency of 6. 00 
cycles per second could then be found from the computer potentiometer 
settings. 



The pitch rate gyro was calibrated by placing it on a constant 
speed turntable and by recording the revolutions per minute and the poten- 
tiometer output. In order to decrease the full scale range, a portion at 
each end of the potentiometer was shorted with metallic paint. The final 
calibration curve is shown on Fig. 9. Full scale deflection represents 
8. 5® per second. 




f 



I 






- -r-idsmaa-- ■ 






f 



I 




r 



1 




12 



The gyro dynamic characteristics were determined by deflect- 
ing the rotor and by recording its return to the equilibrium condition. 
The analog computer was used to determine its damping ratio of .0597 
and its natural frequency of 4. 58 cycles per second. 

The angle of attack vane v/as calibrated in two ways. On the 
ground, a protractor was mounted in such a way that the angular deflec- 
tion of the vane could be read. The angle and potentiometer output were 
recorded and are shown on Fig. 10. An in-flight calibration was made 
by using the telemetering receiver to record the vane position. With 
the aircraft in steady level flight, the aircraft attitude was read with a 
bubble level. By flying at varying speeds, a direct plot of angle of 
attack versus percent telemeter full scale is then available. Fig. 10 
shows both calibrations. The ground calibration indicated telemetering 
full scale was 9. 04 degrees while the inflight calibration gave 7.95 de- 
grees for telemetering full scale. The dynamic characteristics of the 
instrument were not determined due to poor data obtained from this 
instrument. 

In addition to the instrument calibration, it was necessary to 
determine the aircraft weight, c.g. location, and moment of inertia. 
Wheel scales were employed to determine the aircraft total weight. 

The aircraft was then placed on two wing jacks and supported at the tciil 
with £in accurate balance scale. The c. g. location was then determined 
from the location of the support points. These measurements v/ere 
made with the pilot in the aircraft and the wheels raised. The test air- 
craft was equipped with a 55 pound sliding lead weight which was 
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mounted in a hollow tube in the fuselage and which could be moved 17.2 
feet through a crank arrangement in the cockpit. The flight tests were 
made with the lead weight in the aft and in a mid position. The center of 
gravity under these conditions is 29.2 and 32. 5% mac. While measuring 
horizontal distances with the aircraft in the flight attitude, the bomb c.g. 
by 3.4 and 1. 1 inches for the aft and mid c.g. positions respectively. 
Because the bomb c. g. is at a reasonable distance vertically from the 
aircraft c. g. , these distances are dependent upon an accurate attitude 
measurement of the flight conditions at the time of the bomb drop. This 
was accomplished by a flight with the bomb installed using a bubble level. 
The aircraft was then placed in the same attitude while on jacks on the 
ground and the aircraft and bomb c.g. positions were determined. The 
moment of inertia of the eUrcraft v/as determined by oscillating the air- 
craft about the jack points with a spring mounted on the tail. Fig. 11 
shows the configuration employed. The damping under these conditions 
proved to be negligible. The natural frequency was determined by timing 
50 cycles with a stop watch. Knowing the natural frequency and the 
spring constant of the system, the moment of inertia about the jack points 
was determined. Calibration of the springs showed them to be linear for 
the load range required. The moment of inertia was then transferred to 
the aircraft c.g. The moments of inertia about the aircraft c.g. were 
fovind to be 3536 and 3126 slug ft. ^ with the c.g. in the aft and mid 
positions respectively. 
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PROCEDURE 



The procedures employed in this investigation were basically 
very simple. A steady state flight condition was established; the con- 
trols were locked; and the bomb was dropped. The initial airspeed, 
pressure altitude, and temperature were recorded by the pilot. The 
transient response of the aircraft, as sensed by the instruments, was 
transmitted through the telemetering units and recorded on magnetic 
tape. 



In detail, the procediures were carried out as follows. Since 
it was desired to recover the bomb, a water drop was impractical. In 
order to drop over land with reasonable safety, only low altitude drops 
at less than 500 feet were considered practical. In addition, a high alti- 
tude drop would require a major excavation to recover the bomb. 

Since very smooth air was known to be absolutely necessary, 
it was decided to make all drop flights very near sunrise. With gener- 
ally good weather, the air was smooth below 500 feet at approximately 
every other day. This gave an average of about two satisfactory morn- 
ings a week to attempt the drop. Winds of three to five miles per hour 
rendered conditions unsatisfactory for the low altitude test flight. 

The initial intent was to secure good bomb drop data for three 
c. g. positions at a single airspeed. In addition, it was intended to 
record the aircraft response to a small elevator step for at least one 
set of conditions common with a bomb drop. A total of nine bomb drops 
were made. Due to instrumentation and weather problems, no com- 
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pletely satisfactory drop was accomplished. However, two drops were 
made where it was felt that the vinknowns could be removed v/ith reason- 
able confidence. Also, one elevator step was recorded that was felt to 
be usable. The data for these runs are shown in Figs. 14 through 16. 

The pre -flight checks required for a test flight were practically 
non-existent. It was found desirable to check that the aircraft and instru 
mentation were in good condition the afternoon prior to a test flight. The 
weather predictions were checked and arrangements made to gain access 
to all necessary equipment for the drop on the following morning. 

Prior to take-off, radio contact was made with the ground sta- 
tion to ascertain proper functioning of the telemetering installation. 

After take-off, the lead weight was placed in the desired location, £ind 
the aircraft was trimmed for level flight. When proper power settings 
were established the ground station was advised that a run was being 
commenced. The elevator was locked approximately 30 seconds prior 
to the bomb release with a clamp placed on the control column. If air- 
speed and attitude remained constant for the final 30 seconds and no 
gusts were encountered, the bomb was released. It was found that 200 
to 400 feet was the most practical altitude range for bomb release. 
Higher altitudes made bombing accuracy uncertain. Lower altitudes 
were found to lack the necessary smooth air and required excessive 
care in maintaining safe clearances. 

Once take-off was made, the bomb was dropped regardless of 
any instrumentation or telemetering problems. This required that 
streamlining be remounted on the bomb, but it was felt that it was 
xinwise to make a landing with the bomb installed. 
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DISCUSSION AND RESULTS 



Theory 

In studying Equations of Motion of the longitudinal transient 
response of an aircraft, it is conventional to assume that velocity 
remains constant. Also, the aerodynamic forces can be expressed by 
a Taylor expansion about the eqxxilibrium condition retaining only linear 
terms. Under these conditions Ref. 1 develops the equations of motion 
as follows: 

(1) — a + da - dG * fj (t) 

(2) a da + dG - hd^e . {^{t) 

In addition, normal acceleration can be expressed as: 

— 2 ^ n « dG - da 

The particular forcing function in this investigation can lead 
to some confusion in considering initial conditions and in placing the 
equations in non>dimensional form. In this investigation, the aircraft 
transient response was induced by dropping a weight with the controls 
fixed. Clearly, before the v/eight is dropped, the lift is equal to the 
total weight of the aircraft and bomb. Hov/ever, the mass of the air- 
craft during the transient response does not include the bomb weights 
To avoid confusion, it is best to think of the aircraft configuration for 
t < 0 to be that of the basic aircraft (no bomb) with its c. g. deter- 
mined, plus the effects of the bomb acting in relation to the basic air- 
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craft c.g. If the airplane plus bomb configuration is trimmed at t < 0, 
then releasing the bomb is eqtoivalent to applying an up force at the bomb 
c, go Thus, one forcing function in the lift equation is an up force equal 
to the bomb weight, and in the moment equation, we have a positive 
pitching moment eqiial to the bomb weight times the distance between 
the basic aircraft c, go and the bomb Co g. 



(4) 
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bomb weight 



It is pointed out that the initial angle of attack is prescribed 
by as defined above. However, the aerodynamic time {r) and 

the dimensionless moment of inertia (h) are defined by using the air* 
craft mass during the transient response. It is impossible therefore 
to duplicate all conditions using any forcing function other than that 
described previously. 



In addition, it is necessary to examine the expression for 
normal acceleration. This is derived from the expressions developed 
in Refo 1 : 

% • 

(7) F c n W * mV < G - Cr ) 

Z 



If this expression is changed to non* dimensional form, the 



following results: 
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( 8 ) 


T n W 3 


mV (de - da) 


( 9 ) 


n . 


* (dG-da) 


psv^ 


( 10 ) 


n 


« (d©-da) 



The forcing functions in dimensional form are the bomb 
weight in the lift equation; the bomb weight times the distance between 
the aircraft and the bomb center of gravity in the moment equation. 

In non-dimensional form, 'the forcing functions become: 

For t < to fj(t) * f^t ■ 0 

For t > to 
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All time histories obtained in transient responses are in 



real time 
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(13) a + a - 0 s - 
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These expressions are identical to those obtained using an elevator 
step input except for the definitions of the forcing functions. For an 
elevator step input: 



AC 



m a 
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AC 



m 




It is pointed out again that all conditions cannot be matched for the 
two types of forcing functions. In practical termc. initial angle of 
attack and velocity cannot be identical for the two cases at the same 
density altitudeo 



There is an additional point of view which can be taken about 
the initial conditions of the bomb drop. Equations (1) and (2) are 
perturbation equations^ By considering f^ (t) and f^ (t) zero before 
the bomb drop, this establishes a and dG as zero. In a linear 
perturbation problem where the initial and final solutions are both 
steady state conditions, it is arbitrary as to which is considered the 
undisturbed state. It is equally valid, in the aircraft transient re- 
sponse to a bomb drop, to consider the final steady state as the undis- 
turbed condition with a and d0 both zero. The physical motion 
of the aircraft is a steady pitch rate. 
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This can be more clearly demonstrated with the use of equa- 
tions. The variables can be expressed as the sum of a steady state 
and a transient solution. 



a « a a ^ 
as t 



dGsdG + dG. 

ss t 



Under the previous arrangement when t s O; 



s -a* 

o set 



dG = O de 



ss 



d G. 



If the final steady state is^ defined as the xindisturbed condition, the 
equations become : 



“ss = “^8S * ° a * <»t do . de, 



Under these conditions f^ (t) and f^ (t) are zero after the bomb drop 
and are the negatives of the previous forcing fxmctions before the bomb 
drop. That is, the aircraft before the bomb drop has a force and 
moment acting until release. The initial perturbation angle of attack 
and pitch rate are unknown in this case. An analytic solution of the 
equations will show that the steady state solutions in the first case 
are the negatives of the initial condition in the second case. 



The equations for the steady state solutions can be found from 
equations (1) and (2) for the original arrangement. 
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The initial values for d© and a in the second case are repre- 
sented by identical expressions as above. However, the forcing 
fxinctions change signs. 

In the real aircraft problem, there is some disadvantage in 
considering the final state as the imdisturbed condition. The aircraft 
may not reach a complete steady state in the time period where veloc 
ity is essentially unchanged. There is also a mental difficulty in 
thinking of a steady pitch rate as an undisturbed condition in relation 
to level flight. 



l>inear Dependence 

The determination of unique values for stability derivatives 
from flight test data is a problem that arises in any experimental pro 
gram. The stability derivatives which can be determined taniquely is 
controlled by the linear dependence of the equations describing the 
system. This subject is considered in detail in Refs. Z and 3. For 
the longitudinal motions considered in this report, the problem can 
be resolved without great difficulty. 



In equations (13) and (14), it is considered that all flight 

variables and the forcing functions are known and are not zero. 

Under these circumstances, there is no linear dependency between 

the two equations. All stability derivatives and the dimensionless 

moment of inertia (h) can theoretically be determined. The ability 

to determine C and C independently by dropping a 

‘"da ‘"d© 

weight is the basic reason this investigation was undertaken. The 
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use of the elevator as a forcing function does not accomplish this for 
one reason. The tail force introduced into the lift equation is so small 
that its effect on the aircraft motion is negligible. Equation (13) may 
effectively be equated to zero and a relation between the variables a , 
a and 6 results. Therefore, the variables in equation (14) are not 
linearly independent, and it is impossible to find all the stability deriv- 
atives in this equation uniquely. A complete development of what can 
be fouxid is contained in Rel. 2. 

In using a bomb drop as a forcing function, there is an addi- 
tional advantage. The exact nature of the forcing functions are known 
as a force and a moment. With an elevator step, this is not the case. 

It is pointed out that the moment produced by the bomb release does not 
have to be known to solve the problem if h is known. 

The facts stated above can be demonstrated by the use of the 
equations as follows. If equation (13) is used to eliminate 6 from 
equation (14), the following results: 
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If and ^ are Imown and are not zero, all quantities in 

# « # 

parenthesis can be determined with a time history oi a , a and a . 

If h is known, it is not necessary to know the value of ^^m to 
determine these quantities. Inspection of the quantities in parenthesis 
shows that all derivatives are determined uniquely. However, if 
ACl , 

0, as is the case for an elevator input, the derivatives cannot 
be found uniquely even If h and '^m are both known. It is also 
pointed out that ^^m can be identically zero in the case of the bomb 
drop without destroying the uniqueness of the solution. 

Preliminary Investigation 

The previous analysis verified that a solution was possible for 
all longitudinal stability derivatives by dropping a weight. However, it 
was still necessary to verify that the transient response of the aircraft 
would be of sufficient magnitude that the instrumentation available could 
sense the flight variables with the necessary accuracy. It was deter- 
mined that about 500 pounds could be carried externally on the test air- 
craft. This represented approximately 20% of the basic aircraft weight. 
It was estimated that the bomb and aircraft centers of gravity would be 
approximately eight inches apart. Refs. 4 and 5 provided reasonable 
estimates for the longitudinal stability derivatives. With this informa- 
tion, it was possible to estimate the response of the aircraft when the 



bomb was dropped. 
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Ajialog of the Test Aircraft 

It was decided that the use of an anctlog computer was the most 
convenient method to make a preliminary analysis. Equations (13) and 
(14) were set up on the computer and the time history of the aircraft 
simulated flight parameters were placed on a Sanborn recorder. In 
addition, equation (7) was set up to give normal acceleration. See Fig. 
12 for a schematic wiring diagram of the analog computer setup. This 
arrangement allowed variation of the forcing functions and all stability 
derivatives as desired. This arrangemeiit assumed the aircraft was 
perfectly described by the equations and that the instruments would per* 
fectly record the flight variables a , Q and n. 

Solution for Stability Derivatives from Analog Data 

It was found that the statdlity derivatives could be found with 
reasonably good accuracy using the very simple statistical methods dis> 
cussed in Ref. 6. However, it was found that greater accuracy in much 
less time could be achieved. by curve matching two computer solutions. 
One person would set in arbitrary stability derivatives and produce a 
time history of aircraft motions. A second person would then reproduce 
the curves as closely as possible and determine the stability derivatives 

from the potentiometer settings . Excellent results were obtained in this 

C C 

manner for the derivative ^d9 * Good results were found for m^^ 

c c 

while only fair accuracy could be expected for and • 




I 
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It was also discovered that the variation in angle of attack 
would be quite small. This is of cotirse related to the difficulties in 
determining a and a . Improved results could be obtained 
by increasing that is, moving the bomb further from the air- 

craft center of gravity. Unfortunately, the location of the bomb on the 
aircraft was quite closely fixed by the landing gear wheel wells. The 
remaining solution was to shift the aircraft center of gravity. This is 
not nearly as satisfactory aince. in practice, shifting weight ait in the 
aircraft increases its moment of inertia and rapidly decreases . 

Both of these effects tend to reduce the size of the aircraft's response. 
The analog computer demonstrated that moving the aircraft c.g. ait 
increased the response qualitatively about half as much as moving the 
bomb c.g. forward. The character of the changes were distinct how- 
ever. 



In addition, the use of elevator input to create a sizeable 
moment input was investigated. It was found in this case that the re- 

c c 

sponse was effected most strongly by and . It was 

known from the previous theoreticcd discussion that no unique solution 
existed for and • However , the computer 

demonstrated that if a reasonable approximation was available for rate 
derivatives, good results could be obtained for 




as well as 
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c c c 

Using both type inputs, » ”^0! ^ce could be 

c 

reproduced to a close percentage of the initial value. ^da could 

not be reproduced with a high percentage of accuracy, but in terms of 

absolute value. It was well defined. This difficulty with ^da is 

somewhat inherent due to its small size compared to ^m^jg and 

. As Ref. 7 states, the accuracy with which a stability derivative 

can be found is determined by the magnitude of its influence on the motion 

C 

of the aircraft. ’^dO! has a relatively small influence on the air- 
craft's motions and therefore will not be found with high accuracy. 

The variables which were to be measured included normal 
acceleration, pitch rate and angle of attack. Any two of the three will 
provide a complete solution as is shown by equations (1), (2), and (3). 
Under the conditions set into the analog computer for the bomb drop, 
the time histories indicated the following change from the steady state 
in the variables to be measured: 

Normal acceleration: 0 to .4 g 

» 

Pitch rate: 0 to 6® per second 

Angle of attack: -.5*^ to 4-2° 

It was felt that acceleration and pitcli rate could be measured 
in this range with reasonable accuracy. However, it was felt that 
angle of attack could probably not be measured with sufficient accuracy 
to give any useful information. 
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Results of Preliminary lavestigation 

The following facts were indicated by the use of the analog com- 
puter in the preliminary investigation. 

a) The longitodinal stability derivatives could be found uniquely by 
dropping a weight approximately 20% of the aircraft weight. 

b) The weight should be located further from the aircraft center of 

gravity. A distance of 1.00 to 1.75 feet would be satisfactory for the 

Navion. (fk * .176 to. 305) 
c 

Q 

c) The derivative ’^dO would be determined most accurately, 

^^da would be found with only fair percentage accuracy but good abso- 
lute accuracy. The use of an elevator step would assist in determining 
^mQ, and • However, as mentioned previously, all conditions 

of flight cannot be duplicated for the two forcing functions. 

d) It was likely that no reliable angle of attack information would 
be obtained. 

e) The best means of obtaining a solution was probably matcliing 
analog data with the measured response of the aircraft. 

The general technique for finding a set of stability derivatives 
by curve matching was believed to be as follows; 

a) With the best available estimates, create time histories with 
the analog computer of all flight variables which are available from the 
test flight. The step forcing functions are easily created by using the 
reference voltage of the computer. 
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b) Vary each stability derivative separately and make an overlay 
for each derivative showing its effect on the time histories. 

c) Carry out this procedure for each type of forcing function 
considered. 

d) Use the bomb drop time history and create the best match 
possible. 

Cm 

e) Shift to an elevator step. Leave the derivatives d0 and 
^m^jQ fixed, and create the best match possible. 

f) Repeat steps (d) and (e) until the best possible xnatch is 
achieved. 

Some small allowance could be made for the difference in flight 
conditions between the two forcing funetdons. If identical indicated air 
speeds are used, the initial angle of attack for the elevator step will be 
approodmately two degrees less than that for the bomb drop. This con- 
dition is largely relieved by the fact that the angle of attack increases 
three to four degrees during the elevator step while it shows very small 
change during the bomb drop. 

Results Using Flight Test Data 

As previously mentioned in the report, no completely satisfac- 
tory data was obtained. Two bomb drops and one elevator step were 
selected as providing usable information. Figs. 14 through 16 show 
this data. The flight conditions for these runs are recorded in Table 11. 
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When the bomb rack installation was completed, it was found 
that the aircraft and bomb centers of gravity were considerably closer 
together than had been initially estimated. This had a twofold disadvan* 
tage. As previously discussed, it reduces the size of the aircraft 
response, making the Instrumentation requirements more severe. 

Also it makes the determination of the size of the moment forcing func- 
tion more difficult. The horizontal distance between the aircraift and 
bomb centers of gravity depend quite strongly on the aircraft attitude 
when the centers of gravity are very close together horizontally in level 
flight. A one degree error in aircraft attitude created approximately a 
E0% error in the size of the pitching moment when the aircraft c.g. 
was completely ait. The situation was even worse as the aircraft cen- 
ter of gravity moved forward. Since no instrument recorded flight atti- 
tude. it was felt that no accurate moment forcing function could be 
determined. Measurement of aircraft attitude with the bomb installed 
showed a variation too large to be ignored. It was therefore decided to 
regard the momeaat input as an unknown. As brought out in the theory, 
the problem still has a unique solution; however practical application 
showed it was much more difficult. 

The streamlining on the bomb was initially installed to reduce 
the aerodynamic drag. This was thought to be desirable in order to 
estimate the pitching moment more accurately. An estimate of . 5 for 
the drag coefficient of the streamlined bomb was made based on its 
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cross sectional area. This estimate was taken from Ref. 1. The 
moment resulting from the bomb drag was determined to be quite small 
and therefore unimportant since the pitching moment was not accurately 
known anyway. The streamlining was retained throughout this investi- 
gation in an effort to reproduce results. It is not considered to be 
necessary for any future investigation of this nature using low speed air- 
craft. It should be pointed out that it has been assumed that there is no 
aerodynamic effect on the aircraft due to the bomb's presence. 

The matching of the time histories for two bomb drops and one 
elevator step are shown in Figs. 17 through 19. The stability deriva- 
tives obtained through these curves are recorded in Table III. 

In the analog computer solution, the dynamic characteristics 
of the accelerometer and pitch rate gyro were included. The computer 
circuit diagrams for these instruments are shown in Fig. 13. The use 
of filters to remove noise from the test data was considered. In gen- 
eral, it was found that the stray noise could be removed by eye better 
than by filtering when curve matching was employed as a method of 
solution. 



The effect of each stability derivative on the computer time 
histories is shown in Figs. 20 through 24. It can be seen that ^^m, 
and a^^e not clearly defined by the use of the bomb drop. 
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However, and "'da appear to be reasonably well defined. 

c c 

The elevator step improve© the aj>d , but the failure to 

reproduce flight conditions prevents high accuracy. Figs. 25 through 
29 show the effect of each stability derivative on the response to an 
elevator step. 



In spite of the Instrumentation difficulties previously discussed, 
the numerical results of this Investigation are considered to be reason- 
ably good. In the elevator step shown in Fig. 19, there are no ambigui- 
ties in angle of attack or pitch rate. Table IV records the derivatives 
as established in this investigation and those arrived at in Refs. 4 and 5. 

c c c 

The derivatives ^'dG , and ^CL established in this inves- 

tigation fall within the band of previous evaluations. Conditions were 

not reproduced sufficiently to allow a good comparison for . 

Cm 

The indications are that 0i as established in this report is at 

least 10% high. ^O! was determined with the aircraft c. g. located 

at 29 . 2 % and 32,5% mac. This gives a theoretical decrease of . 17 in 

Cm 

the magnitude of . The experimental decrease was .22. 



The calculated and experimental values of ^ for the 

bomb drops agreed within 15%. In both cases, the calculated values 
were greater. See Table 111. 
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Proposed Improvements 

The basic improvements which should be made have already 
been discussed. It is an obvious fact that the instruments employed 
should be reliable and cover the desired range without modification. If 
telemetering is employed, it is a great advantage to use the telemetering 
full scale voltage as the reference voltage for all potentiometer trans- 
ducers. 



It would be of considerable assistance if aircraft attitude and 
elevator position were also measured. While neither is required, they 
both help to determine the quality of the steady state established prior to 
the bomb release. La addition, aircraft attitude could be used to estab- 
lish the moment forcing function. Elevator position would verify that the 
controls remained freed during the bomb release and could be used as a 
forcing function if an elevator input was employed. This would also 
allow the determination of 6 rather than ^ ^m for the elevator 
input. 

The last but most important improvement has already been dis- 
cussed in some detail; the bomb should be moved further from the air- 
craft center of gravity. The importance of this change is shown in 
Figs. 30 through 33 which present analog computer time histories of 
pitch rate and normal acceleration. The moment input here corresponds 
to a distance of 1. 0 feet between the aircraft and bomb centers of gravity. 
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The effect of each stability derivative on tl^iis type of response is shown. 
In these comparisons, the same numerical values were used as those 
employed with the flight test comparisons shown in Figs. 20 through 24. 
With larger distance from the bomb to the aircraft c.g. , the moment 
input could be determined accurately. These curves show that each 
stability derivative has a larger effect on the transient response when a 
larger moment input is used. The relative importance of each derlva- 

Crr. 

tive is not greatly effected; •'however *^0! appears to show the 
greatest improvement. 

There is a characteristic of tliese curves which is of consider- 
able assistance in determining the derivatives by curve matching. 

c c 

There is a short time where ^O! and ”^0! have no effect on the 

transient response. This is represented on the curves by the crossing 

C Ct 

of the variables n and d9 as and are varied. 

This means physically that the aircraft is passing through zero angle of 

attack. Fig. 34 shows that angle of attack initially decreases and then 

increases through zero to a positive steady state value. In the vicinity 

C 

of the minimum, has no effect on the transient response while 

c c 

rrify and Lq, have no effect where a passes through zero. In 
the flight teste, the steady state angle of attack was negative and no 
zero value occurred. 



If the forcing functions are known accurately, approximate 
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values for the stability derivatives will locate the time periods where 

Oi is very small and where the minimum occurs (da = 0), These 

regions may then be used to determine *"d0 and ™da independ- 
- Ct j 

ent of and . 

Ct 

The use of an elevator step to determine a initially is 

advantageous. SiiKe an angle of attack vane can be used with this input, 

C, 

time histories of n aixi a give ^ci directly. For the elevator 
step, equations (7) and (13) show that; 

Cl 

^ a 

V 2 T 

C T 

With known, the remaining derivatives can be determined with 

the analog computer as follows. If a reasonably good match is not avail- 

Ct 

able initially, vary all stability derivatives except O' to achieve a 

Cm C-v, 

fair correlation throughout. Leave ^da fixed and vary "^dO 
Cm 

and ^0! to match the regions where da is very small (the minl- 

Cm 

mum and the steady state of a ). Then vary da to match the 
initial portions of the time histories. 

When a good match is achieved, shift to the elevator step input. 

c c c c 

Leave ^da and ^d0 fixed and adjust ^ a and ^a to 

achieve the best match for all three variables (n» a and dO). Return 

Ct Cm 

to the bomb input and leave a and ™a fixed and achieve the 
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best match possible varying ^da and ^dO . Shift between the 
two forcing functions using these techniques until the best over-all 
match is achieved. 

"While the bomb drop establishes all derivatives uniquely, the 

use of the elevator step will considerably improve the accuracy with 
Cm Ct 

which 0! and 0! can be determined. The accuracy of the 

rate derivatives is improved only slightly. Since an angle of attack 
vane can be used successfully with an elevator step, an extra cross 
check is available. 



36 



CONCLUSIONS 



The proposed method of uniquely determining longitudinal 
stability derivatives for slow speed aircraft is believed to be feasible. 
Due to instrumentation difficulties, positive verification of this is lack- 
ing. However, the numerical results for tlae Navion contained in Table 
m are considered reasonably valid. 



The stability derivatives were determined uniquely from simu- 
lated data generated on tlie analog computer. The simplest and most 
accurate results were obtained by a curve matching technique. 

The use of an elevator step function in addition to a weight 
release v/ill greatly assist in the determination of the longitudinal sta- 
bility derivatives. No consideration was given to any otlier types of 
forcing functions. 



The accuracy with which each stability derivative is determined 
is associated with its influence on the aircraft response. The derivative 

Q 

most accurately determined by the release of a v./eight is ^dO . 



^dc is uniquely found with fair accuracy by a bomb release. 

c c 

and can be found with better accuracy using £in elevator 

^m ^m 

step and previously determined values for dO and da . 
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RECOMMENDATIONS 



As a result of this investigation, the following recommenda- 
tions are proposed: 

1) An investigation of tliis nature should be continued to obtain 
improved experimental results. The configuration and the velocity of 
tlae aircraft should be varied in any additional investigation. 

2) A normal accelerometer with a full scale range of approxi- 
mately one g should be installed in an aircraft for longitudinal 
transient flight testing. 

3) The aircraft should be instrumented to measure elevator 
position and attitude to properly analyse the response from a bomb 
release. 

4) In this type of investigation in a Navion aircreift, a 500 pound 
weight should be at least one foot horizontally from the aircraft center 



of gravity. 
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TABLE I 

Physical Characteristics of the Navion 



A . W ing 

1. Total Wing Area (includes flaps, ailerons 



and 19. 8 ft^ of fuselage) 184.2 ft^ 

2. Span 33.38 ft 

3. Mean Aerodynamic Chord 5.70 ft 

4. Angle of Incidence 

Root +2*^ 

Tip -lo 



5. Twist 

Aerodynamic 

Geometric 

6. Airfoil Section 

Root 

Tip 

7. Aspect Ratio 

8. Taper Ratio 

9. Dihedral 

10. Chord 

Root 

Tip 



B. Vertical Tail 

1. Total Area 

2. Stabilizer Area 

3. Rudder Area 



2° 31' 

3° 

NACA 4415R 
NACA 641 OR 

6.044 

0. 5265 

7. 5° 

7.2 ft 
3.92 ft 

12.93 ft^ 
6. 87 ft^ 
6.05 ft^ 
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TABLE I (Cont.) 



4. Airfoil Section 

Root NACA 0013.2 Mod. 

Tip NACA 0012-64 Mod. 



5. Angle of Stabilizer with fuselage 



2® Nose left 



Horizontal Tail 




1 . 


Total Area 


43.05 ft^ 


2 . 


Stabilizer Area . 


28.95 ft^ 


3. 


Elevator Area 


14.10 ft^ 


4. 


Span 


13. 172 ft^ 


5. 


Mean Aerodynamic Chord 


3 . 34 ft^ 


6 . 


Angle of Incidence 


.40 


7. 


Airfoil Section 


NACA 0012 


8 . 


Aspect Ratio 


4.02 


9. 


Taper Ratio 


.67 


General 




1 . 


Length Overall 


27.25 ft 


2 . 


Tail Length 


16.88 ft 


3. 


Weight (full fuel) 


2430 pounds 


4. 


Engine 


Continental E-185 
185 HP at 2300 RPM, 



29" Hg MoP. at sea level 



5 ? 
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TABLE n 

Test Flight Conditions 





Flight 4 


Flight 9 


Flight 9B 


Date 


23 April 59 


2 May 59 


2 May 59 


Airspeed mph 


94 


96 


94 


Temp 


7°C. 


9°C. 


9°C. 


Pressure 
Altitude (ft) 


350 


150 


450 


c. g. 


29*2% mac 


32. 5% mac 


32 . 5% mac 


Vi (fps) 


137.8 


140. 5 


137.8 


V (fps) 


136.6 


139.4 


137.3 


or 


1.018 


1.0175 


1.006 


W 


2540 


2538 


2525 


o 

O 


.748 


.720 


- 


O 

o 


.611 


. 589 


.607 


A Cl 


. 1360 


. 1312 




I 


3132 


3541 


3541 


h 


.0787 


.0887 


.0887 


T 


1.296 


1.273 


1.298 
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TABLE m 

Stability Derivatives from Flight Tests 





Flight 4 


Flight 9 and 9B 


Date 


23 April 59 


2 May 59 


c.g. location 


29-2% mac 


32 . 5% mac 


C T 


4.98 


4.98 


mQ, 


-.825- 


- .60 


^^dG! 


- .072 


- .073 




- . 167 


- . 166 


A 

m 


.00244 


.00463 (Flight 9. bomb drop) 


Calculated 
for Bomb Drop 


.0029 


.0065 




i 






Comparison of Stability Derivative 
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Fig. 2 Bomb Rack Installation 




Fig. 3 Test Aircraft 
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Fig. 4 Tape Recorder and Telemeter Receiver 




Fig. 5 Analog Computer and Sanborn Recorder 
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Fig. 14 Bomb Drop Flight Data Flight 4 c.g. 29.2% 




Fig. 15 Bomb Drop Flight Data FUght 9 c.g. 32.5% 




Fig. 16 Elevator Step Flight Data Flight 9B c.g. 32.5% 
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